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1 INTRODUCTION 



The Jordan-Brans-Dicke (JBD) |Brans fc Dickdll96ll ) the- 
ory is one of the earliest and most well-motivated alter- 
natives to the Einstein theory of gravitation. The value 
of the gravitational "constant" is set by the inverse of a 
time-dependent classical scalar field with a coupling pa- 
rameter ui. General relativity is recovered in the limit of 
to — > oo. Solar system observations im pose lower bounds 
on (j (|Bertotti et all 120031 : IWilll 120061 ). Generalized JBD 
models are obtained in the low en ergy limit of higher di- 
mens ional theories. Strin g theoretic (iGasperini fc Veneziand 
120071 ) and Kaluza-Klein (|Appelquist et al.ll 19871 ) models af- 
ter compactification of the extra dimensions yield several 
variants of JBD models (or general scalar-tensor models) in 
which the the scalar field coupling to may become dynami- 
cal, and also models with potential terms for the JBD scalar 
field. 

JBD models have been used for tackling several cos- 
mological problems pertaining to different eras of evolu- 
tion of the Friedman-Robertson- W alker (FRW) universe. 
The scenario of extended in flation |Ma thiazhag an fc Johril 
1 19841 ; lLa fc Steinhardtl ll989) was proposed within the con- 
text of JBD cosmology. Generalized JBD models result- 
ing from the compactification of higher dimensional actions 
jMajumdar fc Sethill993 ; IMaiumdar et al.lll993l ; iMaiumdarl 
1 19971 ) enable efficient transition to the post-inflationary ra- 
diation dominated phase (or reheating). The JBD scalar 
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ABSTRACT 

We consider the evolution of primordial black holes in a generalized Jordan-Brans- 
Dicke cosmological model where both the Brans-Dicke scalar field and its coupling to 
gravity are dynamical functions determined from the evolution equations. The evap- 
oration rate for the black holes is different from the case in standard cosmology. We 
show that the accretion of radiation can proceed effectively in the radiation dominated 
era. It follows that the black hole lifetime shortens for low initial mass, but increases 
considerably for larger initial mass, thus providing a mechanism for the survival of 
primordial black holes as candidates of dark matter. We derive a cut-off value for 
the initial black hole mass, below which primordial black holes evaporate out in the 
radiation dominated era, and above which they survive beyond the present era. 
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field has also been incorporated in quintessence scenar- 
ios fo r obtaining the present a cceleration of the uni- 



verse llBertolami fc Martind l200ol ; iBaneriee fc Pavonl [200ll : 
ISen fc Seshadril 120031 ). A more general class of JBD the- 
ory where the coupling parameter u is an arbitrary func- 
tion of the scalar field jBergmannl Il968l ; IWagonerl Il97d ; 
iBergmannl Il97d) has interestin g cons equences on cosmic 
evolution i|Sahoo fc Singh! |2002| . 120031 ) and has also been 
applied in the context of o btaining singularity free cos- 
mology jKalvana Ramal 19971) and relic g ravitational waves 
|Barrow et al.lll993l ; ISingh fc Sahodl2004l ). 

Primordial black holes (PBH s) are pot entially fasci- 
nating cosmo-archeological tools (|Carrll2003t) . It is possi- 
ble for PBHs to impact through their evaporation prod- 
ucts diverse cosmologic al processe s such as baryogene- 
sis and nucleosynthe s is dCarrl 1 19761; IZeldovich et al. I I1977I: 



Mivama & Sate 



syntnes 
d\ 19781; 



19821; iBarrow et al l 



Chcchct kin et al.l. _ _ 
19911 ; IMaiumdar et all Il995l ; iLiddle and Greenl Il998l), the 
cosmic microwave b ackgr ound radiation ( Cha pliriel Il975l ; 
iMacGibbon fc Can _ 1991 ), and the growth of pertur- 



bations as well (|Afshordi et all 120031). PBHs could act 
as seeds for structure formation ( Mack et al.l 120061 ) and 
could also form a significant component of dark mat- 
ter through efficient early accretion in braneworld scenar- 
10s (lMaiumdarll2003l; I Guedens et al.ll20 02; Mai umdar et al.l 
120051 ; IMaiumdar and Mukheried 120051 ). Mechanisms for 
growth of supermassive bla ck holes by PBHs accret ing dark 
energy have been proposed |Bean fc Mag uciio 20 02j) , though 
recent results bv lHarada fc Can (|2005al ) and Harada et al.l 
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(2006) indicate the lack of self-similar growth of black holes 
by accreting quintessence. 

The feasibilit y of black hole s olutions in JBD theory was 
first discussed by iHawkind (1 19721 ). The coexistence of black 
holes with a long range scalar field in cosmology has interest- 
ing consequences such as the possibility of en ergy exchange 
betw een the black holes and the scalar field (|Zloshchastievl 
2005). The black holes themselves could be used to constrain 
the variation of fundamental constants (|MacGibbonl l20(37h . 
Another interesting issue of gravitational me mory of black 
holes in JBD theory has also been studied (|Harada et al.l 
2002). The evolution of PBHs in JBD cosmology has how- 
ever remained unexplored in the literature. Since JBD mod- 
els entail modification of the standard cosmological evolu- 
tion, it is expected that the evolution of PBHs in JBD mod- 
els and their associated consequences on various cosmolog- 
ical processes could be substantially modified compared to 
the case of standard cosmology. 

The purpose of this work is to perform a preliminary 
study of the evolution of PBHs in JBD co smology. We con- 
sider a generalized scalar tens or theory l|Bergmannl 1 19681; 
IWagonerlll970l ; rBergmarmlll970l ) where the coupling param- 
eter u)(4>) is a function of the scalar field <f), and is thus also 
time- dependent. We use a set of power-law solutions for the 
scale factor, the J BD field and the coupling ui obtained in 
the radiation era llSahoo fc Singhll2002l . l2003l ). and investi- 
gate the net energy flux for the PBHs resulting from the 
processes of Hawking evaporation and the acretion of radi- 
ation. We show that accretion of radiation can be effective 
in the radiation dominated era since the growth of PBHs is 
always subdominant to the growth of the cosmological hori- 
zon in this JBD model. In the matter dominated era, the 
dynamics of the PBHs is solely governed by the evaporation 
process whose rate is determined by another set of power- 
law solutions of this JBD mod el in the matter dominated 
era (|Sahoo fc Singhll2002l . [200l ). The lifetime for PBHs are 
obtained in terms of their initial masses. Our results outline 
a viable scenario for the growth and survival of PBHs as 
constituents of dark matter in JBD cosmology. 



2 GENERALIZED BRANS-DICKE MODEL 

We begin with a generalized JBD action given by 

s=^- I d A x^m - ^p-(d^) 2 ] + s mattcr (1) 

I07T J <p 

where the strength of the running gravitational coupling is 
given by the inverse of 0, i.e., at any time t, G — l/4>(t), and 
Smatter corresponds to the action of the relativistic fluid of 
particles (p = p/3) in the flat (fc = 0) FRW early universe 
with scale factor a. The Friedman equation and the equation 
of motion for the JBD field (f> obtained from the above action 
are given respectively by 



p + 3-(p + p) = 



(4) 



a 2 a<j) uj(f 2 p 

n T 7 



60 2 3<j> 



and 



> + 3- 



3p 



ijjd) 



a 2lo + 3 2lo + 3 
with the energy conservation equation given by 



(2) 



(3) 



For a radiation dominated evolution, one has p ~ pn oc 
a~ 4 , and p = p/3. Assuming power law solutions for the 
scale factor a(t), the JBD field <f>(t), and also the coupling 
parameter oj(t) one can obtain (|Sahoo fc Singhll2002l . 120031 ) 
the following time-dependences of these quantities: 



a(t) = di 



4>(t) = & 



_3 

+6 



(•>) 



(6) 



(7) 



where the subscript i indicates the initial values of the vari- 
ables. In order to minimize the departure of the evolution 
of the scale factor from the case of the radiation dominated 
evolution of standard cosmology (i.e., a(t) oc t 1 ^ 2 ) we choose 
\tui\ <C 1. Without loss of generality, we set tfii — £Mpj, 
where M p i = G^ 1 ^ 2 is the Planck mass (corresponding to the 
present value of G) . Note that the solution for the JBD field 
(f> indicates the strengthening of gravity with the increase 
of cosmic time. Such a behaviour for <j> arises not just from 
the specific model ((TJ used here, but is a rather generic fea- 
ture of several generalized JBD or scalar tensor models with 
a potential for the scalar field which r olls down the slope 
leading to the increase of G with time JMaiumdar fc Sethil 
Il992l ; IMaiumdar et al.ll 19931 ; lMaiumdar|]i997l ). 



3 EVOLUTION OF PRIMORDIAL BLACK 
HOLES 

Primo rdial black holes could form due to various mecha- 
nisms (|Carrll 19851) some of which have been studied in detail 
I" 



in the literature (fZeldovich fc Novikov 1966; Hawking il97ll: 
Carr fc Hawking] 1 19741 ; lHawking et al.|[l982l ; iKodama et al 



19821 ). The gravitational collapse of perturbations within the 



cosmological horizon could lead to sub-horizon sized black 
hole formation after Jeans crossing. Such a mechanism may 
not be very effective however, in JBD cosmology where the 
Jeans length itself could be significantly increased due to the 
weakening of gravity in certain models by a large value of 
the JBD field in the early universe. Nevertheless, a strong 
(first order) inflationary phase trans ition in extended infla- 
tion models driven by the JBD field (lLa fc SteinhardtlllQSL 
Maiumdar fc Sethilll992l ; IMaiumdar et al.lll993l ; iMaiumda: 



19971 ) naturally leads to unnucleated and trapped false vac- 



uum regions and topological defects such as domain walls 
and wormholes that could easily collapse into black holes, 
and could also create super-horizon density perturbations. 
The formation of super-horizon scale PBHs i n the expanding 
FRW b ackground has been studied recently l|Harada fc Carr] 
l2005bl ). In the present analysis we will not take recourse to 
any particular formation mechanism for PBHs, but rather 
study their cosmological evolution in JBD theory, assuming 
that there exist PBHs in such scenarios. 

We now consider the evolution of PBHs in the cosmo- 
logical background governed by the above solutions (O,© 
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and 0. We assume that the PBH density is low enough to 
ensure radiation domination. For a PBH immersed in the ra- 
diation field, the accretion of radiation leads to the increase 
of its mass with the rate given by 



M acc = 4irfr BH p R 



(8) 



where tbh = 2M/cj> is the black hole radius, and / is 
the accretion efficiency. Using the solution for <f> given by 
Eq.© with the assumption of \u>i\ <C 1, and using pR = 
[(3M^)/(327rf 2 )], one obtains 



M a cc = 



BM 
t 



(9) 



where B = (3/)/(2£ 2 M 2 ; ^)- Eq.© is integrated to yield 

M 1 

M 



(10) 



1 - BMM(t/ti) 

with Mi being the initial mass of the PBH at time t%. Since 
the logarithmic growth rate for a PBH given by Eq.fTUJ) 
is subdominant to the linear growth of the horizon mass 
Mh ~ t (since a ~ t 1//2 ), once a PBH is formed, it is indeed 
possible for it to grow in size by accreting the radiation 
energy within its cosmological horizon. 

For a complete picture of PBH evolution, one needs also 
to consider the Hawking evaporation process, whose rate is 
given by 



Mcva 



-4:Tvgr 2 BH aT 4 



(11) 



where T = (j>/(8irM) is the Hawking temperature, a the 
Stefan-Boltzmann constant and g is the effective number of 
degrees of freedom of the particles emitted by the black hole. 
The solution for <f> given by Eq.© leads to 

A 



M, 



MH 



(12) 



where A = (ga^ 2 M 4 t ti)/ (2567T 3 ). The complete evolution 
for the PBH is thus described by the combination of Eqs.([9} 
and d): 

It is apparent from Eq. (|13[l that for PBHs with initial mass 
Mi « (A/B) 1 / 4 , the rate of evaporation exceeds that of 
accretion. For such a case, accretion soon becomes negligi- 
ble and black holes lose energy at a rate given effectively by 
Eq. (|12[) . Note that though the rate of evaporation decreases 
with time, it is still higher than the corresponding rate in 
standard cosmology. This is because the Hawking temper- 
ature T = 0/(8ttM) is larger for JBD PBHs for large <j>. 
Hence, a PBH with initial mass Mi < (A/B) 1 ^ 4 evaporates 
out much quicker, with a lifetime t evap given by 

3 



t, 



exp 



1 
3l 



Mi 
M v i 



(14) 



However, PBHs with initial mass Mi > (A/B) 1 ^ 4 expe- 
rience monotonic growth with accretion dominating over 
evaporation throughout the period of validity of Eo. (ll3l) . 
i.e., throughout the radiation dominated era. Eq. (|13[) can 
be integrated exactly and leads to the following mass-time 
relationship for the PBHs: 



ln(- 



1 



2BC 1 / 4 



tan 



M 

cm 



tan 



+ 



(M - C 1/4 )(M, + C 1/4 ) 
4BC 1 / 4 (Mi - C" 1 /4)(M + C 1 / 4 ) 1 



-ln| 



(15) 



with C = A/B. Accretion of radiation can proceed ef- 
fectively till the universe stays radiation dominated, i.e., 
up to the era of matter radiation equality t eq . This re- 
sult is qualitatively different from the widely accepted pic- 
ture in the standard cosmological evolution where accre- 
tion of radiation in the rad iation dominated era seems 
to be ineffective (|Card 120031 '). The domination of accre- 
tion over evaporation is observed also in other modi- 
fied gravity theor i es, such as in t he br a neworld scenario 
ilMajumdarl 120031; iGuedens et al.1 |2002| ; iMaiumdar et all 



numd an l/uud; uucdcns et al.l izuua IMaiumaar et au 
3 Maiumdar and Mukheriedl2005l )~The maximum mass 



12005 

achieved by a PBH of initial mass Mj is given by 

R^max 1 



Mi 



BMi\n(t eq /U 



(16) 



In the matter dominated era, pu ~ a 3 and p — 0. A 
set of solutions for t he JBD cosmological equat ions ([2]), ([3jl 
and © is given by (|Sahoo fc Singh 2002, 20o3) 



a(t) = a(t eq ) 



<P(t) = 4>{tc 



u(t) = -w(tesr) ( -r- 



(17) 



(18) 



(19) 



Matching these solutions with those obtained in the radi- 
ation dominated era given by Eqs.([SJ,© and 0, one has 
<t>{teq) = CM^iU/teq) 1 ^ and tj (teq) = uJi(U q /U) 1/2 . Substi- 
tuting the solution for the JBD field <f> in Eo. (|lip leads to 
the following evaporation law for PBHs in the matter dom- 
inated era: 



M = -A 



MH 8 / 3 



(20) 



with A = AtiM^t% 3 . The mass of the PBHs evolve as 



M(t) = Mmaz 



3.4 



-5/3 



t 



-5/3 



1/3 



(21) 



in the matter dominated era. The lifetime for the PBHs t eV a P 
is given by 



tevap — te 



5 M„ 



-3/5 



(22) 



which upon (using Eq. (|16[l ) reduces to 



tevap — te 



Mf 



1 - BMMteq/U)) (AUM*,) 



-3/5 



(23) 



Since the availability of background radiation diminishes 
substantially in the matter dominated era, one can safely 
neglect any accretion term corresponding to the accretion 
of radiation by the PBHs. However, other forms of accre- 
tion such as of the en ergy of a quintessence fi eld could, in 
principle, be effective ()Bean fc Maguei"ioll2002h . that we do 
not consider in the present analysis. The rate of evapora- 
tion is again faster, as in the radiation dominated era, than 
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in the case of standard cosmology (where t evap ~ M^ ax ). 
However, the rate decreases with time as is evident from the 
r.h.s of Eg. pop whose counterpart in the standard theory is 
independent of time. 

The overall lifetime of a JBD PBH has an interesting 
comparison to a PBH lifetime in standard cosmology. A 
JBD PBH with initial mass Mi < (A/B) 1/4 , evaporates out 
quicker than a PBH in standard cosmology, being unable to 
accrete in the radiation dominated era. However, when the 
initial mass exceeds (A/B) 1 ^ 4 , accretion proceeds effectively 
and dominates over evaporation in the radiation dominated 
era increasing the PBH mass at the time of matter-radiation 
equality. In the matter dominated era, though the evapora- 
tion rate of JBD PBH is faster, it starts to evaporate out 
much later than in standard cosmology (without accretion) . 
Thus, the lifetime for a PBH with Mi > (A/B) 1/4 is en- 
hanced in the JBD scenario compared to the case of stan- 
dard cosmology. Note that in the present analysis we have 
not considered any possible back reaction of the PBHs on the 
local background value of the JBD field resulting from the 
local change of energy density p due to the PBHs. The form 
of the JBD field equation © is such that the p-dependence 
drops out in the radiation dominated era. Further, in the 
late matter-dominated era the value of uj is large enough 
(in order to satisfy present observational bounds) such that 
the p-dependent term again becomes negligible in Eq. © . So 
any back reaction could be effective only in the early matter 
dominated epoch during which a local rise in energy density 
could result in a slower rate of decrease for the <f> field near 
the PBHs than what is obtained in Eq. (|18[) . The resultant 
change in the PBH evaporation law (|20[) could introduce 
a minor correction to the PBH lifetime (|23p that we have 
neglected here. 

It would be interesting to observe particular examples 
of evolving PBHs in JBD cosmology. Note that the accretion 
and evaporation rates depend on the initial value <f) of the 
JBD field through the quantities A and B defined earlier. 
Here we fix the value of <f>i or £ by setting the present value 
of <j}(t n0 w) = Mp t (corresponding to the present strength 
of gravity) and evolving (j> given by Eqs. (|18p and © corre- 
sponding to the matter and radiation dominated eras, re- 
spectively, backwards in time. We thus obtain the value of 
4>{Pew) = 10 18 Mpj, where tEW denotes the electroweak era 
since which the JBD evolution equations ([3]) are assumed 
to be valid. From Eq. (fl4)l it follows that for a PBH with 
Mi < (A/B) 1 ' 4 to survive the radiation dominated era, 
Mi > 10 13 M pi . A PBH which forms at time t E w with an 
initial mass of the order of the cosmological horizon mass 
at tEW satisfies the condition Mi > (A/B) 1 ^ 4 , and accretes 
radiation. Though the actual growth of mass turns out to be 
negligible for such a PBH, accretion does play an important 
role enabling it to survive for much longer (t evap > t now ). 

Another interesting consequence of Eq. ([23} is that for a 
PBH to evaporate during the present era, i.e., t evap = t now 
one needs to have M max « 1Q 10 M p i. Since accretion can- 
not cause significant growth, one requires Mi ~ M max w 
10 10 A/ P ;. However, a PBH with such a small initial mass 
cannot even survive the radiation dominated era. Only those 
PBHs which have Mi > 10 13 M P ; enter the matter domi- 
nated era and these also survive up to the present era, i.e., 
t e vap » tnow- Thus there are no end products of black 
hole evaporation to contend with in JBD cosmology during 



the matter dominated era and beyond up to present times. 
All PBHs with initial mass Mi < 10 13 M P ; evaporate out 
in the radiation dominated era itself. On the other hand, 
PBHs with larger mass can survive as intermediate mass 
astrophysical black holes contributing to the dark matter 
density. 



4 CONCLUSIONS 

To summarize, in this letter we have studied the evolution of 
PBHs in JBD cosmology and obtained a cut-off value for the 
initial mass which decides whether the PBHs could survive 
up to today. An interesting consequence of our analysis is 
that no PBHs evaporate during the matter dominated era, 
and hence there is no distortion of the cosmic microwave 
background radiation spectrum due to evaporating PBHs. 
Though our results are in the context of a particular JBD 
model, the sort of evolution obtained for the JBD field </> is 
generic to other JBD models as well, in particular those fol- 
lowing from higher dimensional theories |Maiumdar fc Sethil 
1 19921 ; iMaiumdar et all 1 19931 ; lMaiumdarlll997h . In order to 
satisfy solar system constraints (Will 2006) such as the one 
recently obtained through the frequency shift of radio pho- 
tons to and from the Cassini spacecraft as they passed near 
the sun (|Bertotti et al.ll2003l ). one requires \uj(t now )\ > 10 4 , 



which using Eqs. ([19} and ([?}■ leads to uji ~ 10~ 14 , thus 
strongly validating our approximation of <C 1 used to 
derive the black hole evolution equation in the radiation 
dominated era. With the possibility of future observations 
of act ually finding or ruling out JBD models by determin- 
ing LU (|Acauaviva fc Verdd [20071 ), it would indeed be excit- 
ing to re-work the standard observationa l constraints |Carrl 
120031 ; iLiddle and Greenlll998l ; ICarrlll985l ) on the density of 
PBHs at several cosmological eras in the JBD scenario. In 
particular, the faster evaporation of smaller PBHs in the 
radiation dominated era could impose tighter constraints 
on the initial mass spec trum from nucleosynthesis bounds 
jLiddle and Greenlll998h . Investigation of these issues calls 
for a more comprehensive analysis by considering a popula- 
tion of black holes in the early universe evolving according 
to the JBD dynamics. 
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